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ABSTRACT 

We present photometry and spectra of the type IIP supernova 1999em in NGC 1637 
from several days after the outburst till day 642. The radioactive tail of the recovered 
' bolometric light curve of SN 1999em indicates the amount of the ejected 56 Ni to be 

| w 0.02 Mq. The Ha and He I 10830 A lines at the nebular epoch show that the 

^5 distribution of the bulk of 56 Ni can be represented approximately by a sphere of 56 Ni 

'-^ | with a velocity of 1500 km s -1 , which is shifted towards the far hemisphere by about 

400 km s" 1 . The fine structure of the Ha at the photospheric epoch reminiscent of the 
Q ■ "Bochum event" in SN 1987A is analysed in terms of two plausible models: bi-polar 

Jjj | 56 Ni jets and non-monotonic behaviour of the Ha optical depth combined with the 

c/3 . one-sided 56 Ni ejection. The late time spectra show a dramatic transformation of the 

[O I] 6300 A line profile between days 465 and 510, which we interpret as an effect 
of dust condensation during this period. Late time photometry supports the dust 
formation scenario after day 465. The [O I] line profile suggests that the dust occupies 
a sphere with velocity s=s 800 km s _1 and optical depth 3> 10. The latter exceeds the 
optical depth of the dusty zone in SN 1987A by more than 10 times. Use is made 
of the Expanding Photosphere Method to estimate the distance and the explosion 
time, D w 7.83 Mpc and to ~ 1999 October 24.5 UT, in accord with observational 
constraints on the explosion time and with other results of detailed studies of the 
method (Hamuy et al. 2001; Leonard et al. 2002). The plateau brightness and duration 
combined with the expansion velocity suggest a presupernova radius of 120 — 150 i?©, 
ejecta mass of 10 — 11 M© and explosion energy of (0.5 — 1) x 10 51 erg. The ejecta 
mass combined with the neutron star and a conservative assumption about mass loss 
implies the main sequence progenitor of M ms « 12 — 14 M©. The derived mass range is 
in agreement with the upper limit to the mass found using pre-supernova field images 
by Smartt et al. (2001). From the [O I] 6300,6364 A doublet luminosity we infer the 
oxygen mass to be a factor four lower than in SN 1987A which is consistent with 
the estimated SN 1999em progenitor mass according to nucleosynthesis and stellar 
evolution theory. We note a "second-plateau" behaviour of the light curve after the 
main plateau at the beginning of the radioactive tail. This feature seems to be common 
to SNe IIP with low 56 Ni mass. 

Key words: supernovae: general, line formation, dust formation - supernovae: indi- 
vidual: SN 1999em, SN 1987A - abundances, stars: evolution. 
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1 INTRODUCTION 

Type II (and Ib/Ic) SNe are generally associated with re- 
gions of recent star formation in spiral galaxies, suggest- 
ing that they represent the final episode in the life of mas- 
sive stars (M > 8M Q ) which explode owing to core collapse 
(Filippenko 2001). The study of core collapse events is im- 
portant for understanding the range of progenitor masses 
which produce them, the consequent nucleosynthesis for its 
effect on galactic chemical evolution and the explosion en- 
ergy which remains an ill-determined quantity for the vast 
majority of SNe, and which is also relevant for gas dynami- 
cal processes and ejection of material from galaxies. Detailed 
photometric and spectroscopic observations of SNe II on a 
long time scale are still rare, especially for SNe IIP (plateau). 
Meanwhile from the recent experience with SN 1987A we 
know, how valuable can be extended sets of photometric and 
spectroscopic data for understanding what really happens to 
the massive (M > 8Mq) star when its central iron core col- 
lapses. In this way one obtained information on the amount 
of ejected Ni and oxygen, and their mixing and dumpi- 
ness, possible asymmetry of 56 Ni ejecta, "Bochum event", 
dust formation. The case of SN 1987A also provided us with 
the possibility for testing and revising the theory of stellar 
evolution, SN hydrodynamical models and models of spectra 
formation. 

Supernova 1999em, discovered by the Lick observatory 
Supernova Search on Oct 29 UT at an unfiltered magnitude 
~ 13.5 mag in the nearby galaxy NGC 1637 (D = 7.8 Mpc; 
IAUC 7294), has become another well observed type IIP 
event. Detected very soon after the explosion and followed 
for more that 600 days this event gives another boost to 
studies of SNe II. SN 1999em is the first type IIP detected 
at both radio and X-ray wavelengths at early time (Poo- 
ley et al. 2001). However it is the least radio luminous and 
among the least X-ray luminous SNe (Chandra X-ray and 
radio NRAO observations; IAUC 7318, 7336 ). This early 
and weak radio emission is consistent with a low mass loss 
shortly before the explosion for SNe IIP (Baron et al. 2000) . 
The analysis of the spectra at the early photospheric epoch 
already emphasized the problem of He abundance in the hy- 
drogen envelope of SN IIP progenitors (Baron et al. 2000) 
and permitted one to check and upgrade the method of the 
expanding photosphere (EPM) for the distance determina- 
tion (Hamuy et al. 2001; Leonard et al. 2002). The spec- 
tropolarimetry on the other hand has been studied at 5 dif- 
ferent epochs (till day 163 after discovery; Leonard et al. 
2000). The authors estimated the broadband polarization 
at day 7 to be ~ 0.2 %, rising to ~ 0.5 % on day 161. The 
low polarization found for SN 1999em and hence the small 
implied asphericity especially at early photospheric phase 
is another encouraging reason for the validity of the EPM 
for this class of object (SNe IIP). The progenitor nature of 
SN 1999em was discussed by Smartt et al. (2002), who used 
pre-explosion CFHT images to derive bolometric luminos- 
ity limits and thus constraints on the mass of the progenitor 
star of SN 1999em. They concluded that the main-sequence 
mass should be < 12 ± 1M Q . 

Here we present photometry and spectroscopy of 
SN 1999em from the phase of several days after the explo- 
sion till day 642. In what follows we describe the photometry 
and spectral evolution and provide estimates of the 56 Ni and 
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Figure 1. SN 1999cm in NGC 1637. Image taken at ESO on Mar 
27. 2000, about 150 days after discovery (in the R band). The SN 
position is indicated, also shown are the reference sequence stars. 
North is up, East is to the left. 



oxygen mass in SN 1999em; we will analyse fine structures 
of the line profiles at the photospheric and nebular epoch 
in an attempt to determine the effects of the 56 Ni mixing 
and asymmetry. The last spectra are used as diagnostics of 
possible dust formation in the ejecta. The photometry and 
spectroscopy at the plateau phase will be used to recover 
global parameters of SN (presupernova radius, mass and en- 
ergy of SN). We then discuss the main sequence mass of the 
progenitor and some important correlations. 



2 PHOTOMETRIC EVOLUTION AND 
BOLOMETRIC LIGHT CURVE 

Figure 1 shows the position of SN 1999em in its host 
galaxy NGC 1637. It is located at R.A=4:41:27.13, a Decl=- 
02:51:45.2, about 15" west and 17" south of the nucleus of its 
host galaxy NGC 1637. NGC 1637 is an SABc late type spi- 
ral galaxy having a heliocentric recession velocity of Vh c i ~ 
717 km s _1 . According to Sohn & Davidge (1998) it has a 
distance of 7.8t£;g Mpc. A KAIT* image of the field taken 
on Oct 20.45 showed nothing at the position of SN 1999cm, 
limiting the magnitude to about 19.0, suggesting that it was 
discovered shortly after the explosion. 



2.1 Interstellar extinction 

Galactic extinction towards NGC 1637 is known to be A^? 1 = 
0.134 ( Schlegel et al. 1998) corresponding to a colour excess 
E(B — V)=0.04 (adopting the standard reddening laws of 
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Figure 2. UBVRI Light curves of SN 1999em. The light curves have been shifted by the reported amounts. Also shown is the slope of 
56 Co to 56 Fe decay. Reddening corrections have not been applied. 



Cardelli et al. 1989). However estimating the extinction due 
to the host galaxy is more complicated. 

Recently Baron et al. (2000) have provided constraints 
on the total colour excess using theoretical modelling of two 
early time spectra, namely their 29 Oct & 4-5 Nov spectra. 
They found that E(B-V) = 0.05-0.10 and E(B-V) < 0.15. 
A value of E(B - V)=0.05, as derived from their SYNOW 
fit to their Oct 29 spectrum, implies an amount of E(B — 
V)=0.01 due to the host galaxy corresponding to A\? Bt = 
0.031 and y4g ost = 0.041. Adopting, for example, these values 
one finds M^n* < —15.8, which places SN 1999em at the 
lower extremity of the regular type Il-plateau SNe in the 
classification of Patat et al. (1994). 

On the other hand we note the presence of a narrow 
Na I D interstellar absorption line in our first spectra at the 
redshift of NGC 1637 with an equivalent width (EW) ~ 1.3 
A, and therefore we would expect significant additional red- 
dening in the host galaxy. In fact some correlations in the 
literature relating the equivalent width of Na I D to the red- 
dening within the host galaxy indicate values of the order 
E(B - V) ~ 0.32 (Barbon et al. 1990) and E(B - V) ~ 0.17 
(Benetti et al., unpublished correlation), corresponding to 
0.99 and A\° st = 0.52, respectively. If, for instance, 

16.48, which 



\ host 



Ay = u.yy fclllU I\y 

we assume this last value, we obtain M,^ ax < 



implies that SN 1999em is a regular SN IIP but with sig- 
nificant reddening. Uncertainties in the calibration of ab- 
sorption line strengths prevent a confident determination of 
reddening by this method. 

To overcome this disagreement between theoretical 
modelling and empirical calibration, Hamuy et al. (2001) 
have tested the sensitivity of the Expanding Photosphere 
Method (EPM) to reddening caused by the host galaxy in 
deriving the distance to the supernova (using different filter 
subsets, their Figure 12). They conclude a consistent value 
for the extinction in the host galaxy is Ay st = 0.18 and 
thus a total extinction Ay* = 0.31 and Ag* = 0.41, which 
we adopt as the best estimate for our present analysis. 



2.2 Light curves 

The UBVRI photometry data are reported in Table 1 to- 
gether with the different telescopes and instruments with 
which the observations were obtained. Our photometry 
started at ~9 d and extended to ~508 d (after explosion 
time). The reduction was carried out within the IRAF envi- 
ronment, applying bias, overscan and flat-field corrections. 
The contamination of the galaxy has been removed using 
PSF substraction. The photometric calibration of the SN 
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Table 1. SN 1999cm Photometric Observations 



Date 
(UT) 


JD 
2400000+ 


U(a) 


B(a) 


V(a) 


K(a) 


1(a) 


Instrument 


28/10/99 


51479.51 






13.69 .06 









03/11/99 


51485.62 


13.02 .03 


13.85 .03 


13.84 .03 


13.67 .03 


13.63 .03 


4 


03/11/99 


51485.81 




13.89 .03 


13.84 .03 


13.65 .03 


13.61 .03 


2 


04/11/99 


51486.62 


13.04 .03 


13.90 .03 


13.91 .03 


13.71 .03 


13.59 .03 


4 


04/11/99 


51486.73 




13.93 .03 


13.85 .03 


13.68 .03 


13.61 .06 


2 


05/11/99 


51487.57 


13.07 .12 


13.90 .06 


13.98 .03 


13.83 .09 


13.55 .12 


4 


05/12/99 


51517.80 


15.92 


14.93 


13.98 


13.64 


13.40 


1 


08/12/99 


51520.56 






13.92 .12 






5 


11/12/99 


51523.54 


16.01 .03 


15.01 .03 


13.95 .03 


13.50 .03 


13.30 .03 


4 


13/12/99 


51525.51 


16.04 .03 


15.06 .03 


13.98 .03 


13.52 .03 


13.27 .03 


4 


15/12/99 


51527.57 


16.15 .03 


15.07 .03 


13.99 .03 


13.50 .03 


13.27 .03 


4 


18/12/99 


51530.73 


16.41 .27 


15.11 .09 


14.02 .03 


13.55 .03 


13.29 .03 


1 


25/12/99 


51538.59 


16.64 .06 


15.20 


14.04 


13.54 


13.26 


4 


27/12/99 


51540.50 


16.68 .03 


15.25 


14.05 


13.51 


13.28 


4 


28/12/99 


51541.71 


16.82 .09 


15.28 .06 


14.05 .03 


13.61 .03 


13.24 .03 


2 


05/01/00 


51548.51 


16.87 .03 


15.33 .03 


14.05 .03 


13.55 .03 


13.29 .03 


4 


O"! / r\~i 1 c\r\ 

31/01/00 


51575.44 




15.81 .12 


14.35 .09 


13.80 .03 




5 


09/02/00 


51584.32 




16.22 .39 


14.66 .09 


14.05 .06 


13.64 .06 


5 


16/02/00 


51590.53 


18.55 


16.29 


14.76 


14.11 


13.77 


2 


13/03/00 


51616.50 






16.46 .06 


15.57 .06 




2 


08/04/00 


51642.52 




18.17 


16.66 


15.76 




1 


31/08/00 


51787.71 






18.04 .06 


17.16 .03 




3 


31/08/00 


51787.74 


20.51 .21 


19.13 .06 






16.66 .06 


3 


01/09/00 


51788.90 


20.66 1.05 






17.18 .06 




1 


08/09/00 


51796.60 






18.20 .15 


17.37 .09 




5 


17/11/00 


51866.77 


20.94 .72 


19.81 .54 


18.82 .24 


18.05 .12 


17.62 .24 


2 


01/02/01 


51941.65 


20.85 .12 


20.24 .18 


19.54 .12 


18.86 .12 


18.39 .09 


1 


16/03/01 


51985.53 




20.69 .27 


20.05 .18 


19.47 .24 


18.91 .21 


2 



= WHT(Imagc provided by Smartt.); 1 = ESO 3.6m+EFOSC2; 2 = Danish 1.54m+DFOSC 

3 = TNG+DOLORES; 4 = TNG+OIG; 5 = Asiago 1.82m+AFOSC 
* The reported errors are la. 



was made relative to the local sequence of stars (see Figure 
1) in the field of the host galaxy NGC 1637 and calibrated 
using observations of Landolt standard stars (Landolt 1992), 
obtained when the nights were photometric. The SN mag- 
nitudes have been measured using a PSF-fitting technique. 

The UBVRI light curves of SN 1999em are illustrated 
in Figure 2. The Figure includes also data from Hamuy et 
al. (2001) as well some data from circulars and VSNET 
service^. We have also plotted the decay slope of 56 Co : 0.98 
mag (100 d) _1 , which corresponds to the late-phase decline 
rate of most normal SNe II L-P, especially for the V-band 
(Turatto et al. 1990 and Patat et al. 1994). 

At early times, especially in the V band, one notes that 
SN 1999em rises to a brief peak, then drops by about 0.6 
mag in ~ 3.5 days, then rises again to a brief second max- 
imum followed by a settling onto the plateau phase. Here 
the expansion of the photosphere and the cooling balance 
each other to produce a phase of almost constant magni- 
tude, m v ~ 14 and M v ~ —15.76. The plateau phase lasts 
about 80 days, after which the light curve displays a steep 
decline (~ 2 mag in ~ 40 days for the V-band), signaling the 
onset of the nebular phase and the start of the exponential 
decline. Similar early short duration rapid change was seen 

t http://www.kusastro.kyoto-u.ac.jp/vsnet/ 



in SN IIP 1988A (Ruiz-Lapuente et al. 1990) as well as for 
SN 1993J (Barbon et al. 1995), and can be interpreted to be 
a consequence of low mass loss in the immediate presuper- 
nova phase (see discussion of the bolometric light curve). 

After the steep decline from the plateau phase and start- 
ing around 130 days after explosion the shapes of the light 
curves evolve with time: firstly there is a clear flattening 
from day 130 lasting about 50 days, especially at blue wave- 
lengths, which is then followed by the exponential decline. 
Similar flattening is also seen in the data for the faint SN 
IIP 1997D (Benetti et al. 2001) and is present in the light 
curves of SN IIP 1991G (Blanton et al. 1995). 

A linear fit to the tail (from ~180 to ~ 510) gives 
the following decline rates (in mag (100d) _1 ): jb ~ 0.66, 
7v ~ 0.97, 7h ~ 1.08 and 7/ ~ 1.07 which are close to the 
values found for the typical SN IIP 1969L (Turatto et al. 
1990). Especially for the V-light curve this indicates that 
radioactive decay of 56 Co with consequent trapping of 7- 
rays is the main source of energy powering the light curve 
at late times, at least until ~ 510 days after explosion. 

We note also that the late U light curve remains almost 
flat, similar in behaviour to that reported for SN 1987A until 
~ 400 days after explosion (Suntzeff et al. 1988). 
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Figure 3. B - V colour of SN 1999em and SN 1987A. Bottom 
panel shows the difference A(B — V) = (B — V)99 em — (B — V)s7A ■ 
Both SNe have been corrected for reddening. For SN 1987A we 
use Ay = 0.6 (Filippcnko 1988). Filled circles represent data from 
Hamuy et al. 2001. 



Figure 4. Comparison of the V absolute light curve of SN 1999em 
with those of SN 1987A and SN 1997D. Distances moduli of fi = 
29.46, fi = 18.49 and /i = 30.64 are assumed for the three SNe, 
respectively. Dashed line shows the luminosity dependence due to 
56 Co decay for arbitrary mass of 56 Ni. 



2.3 Colour evolution 

In Figure 3 we show the (B — V) colour curve of SN 1999em, 
together with that of SN 1987A for comparison. The explo- 
sion date of SN 1987A is well established to be February 
23.316 from the Kamiokande and 1MB neutrino detections 
(Hirata et al. 1987; Bionta et al. 1987), while for SN 1999em 
the explosion date is assumed to be October 24.5 (~ JD 
2451476) derived from the EPM, to be discussed later. 

The (B — V) colour varies at different rates accord- 
ing to phase. During the first 40 days, it exhibits a steep 
and rapid decline from blue (high temperatures) to red (low 
temperatures) as the supernova envelope expands and cools. 
Subsequently the (B — V) colour varies more slowly as the 
rate of cooling decreases. At ~130 d the (B — V) colour 
reaches a value around 1.7, then it turns blue again (rate 
~ —0.33 mag (100 d) _1 ), as the light curve settles into the 
exponential tail. 

Note the general similarity of (B— V) colour evolution of 
SN 1999em with that of SN 1987A especially at late phases, 
where the A(B — V) approaches zero (the bottom box in 
Figure 3). This fact gives us confidence that comparison of 
absolute V-light curves should be reliable for the estima- 
tion of 56 Ni mass, released in the explosion of SN 1999em, 
because the bolometric light-curves are liable to be similar 
also. 

There is little doubt that the difference of the colour be- 
haviour of both supernovae at the early epoch (faster cooling 
of SN 1987A) is related to the difference in presupernovae 
radii. The smaller radius of SN 1987A presupernova leads 
to a faster adiabatic cooling of the radiation trapped in the 
expanding envelope. The details of the colour behaviour of 



SN 1999em should be reproduced in future hydrodynamical 
modelling. 



2.4 The bolometric light curve and Jb Ni mass 

In Figure 4, we present the absolute light curve of SN 
1999em in the V band. The light curve of the double-peaked 
SN 1987A and the peculiar type IIP SN 1997D are also dis- 
played for comparison. These plots highlight the different 
behaviour of the photometric evolution of the three SNe, 
from the explosion until about 550 days. They provide con- 
straints on the explosion energy of SN 1999em as well as on 
the radioactive 56 Ni mass ejected. All the radioactive tails 
of the three SNe, especially in the V band, follow the 56 Co 
decay slope suggesting that it is the main energy source pow- 
ering the late phases of the light curves. The fainter object 
SN 1997D, My™ > -14.65 (Turatto et al. 1998), ejected 
an extremely small amount of radioactive 56 Ni, about 0.002 
Mq , with a low explosion energy of ~ 10 50 ergs derived from 
modelling the spectra (Chugai & Utrobin 2000), while the 
mass of 56 Ni ejected by the well studied SN 1987A and the 
explosion energy are ~0.075 Mq and ~ 1.3 x 10 51 ergs, re- 
spectively (Danziger et al. 1988; Woosley et al. 1989). We 
therefore conclude that SN 1999em produced a mass of 56 Ni 
intermediate between 0.002 and 0.075 Mq. To obtain a more 
accurate estimate of the 56 Ni mass we constructed the "bolo- 
metric" light curve of SN 1999em by integrating the flux in 
the UBVRI bands. The resulting light curve is shown in 
Figure 5 together with that of SN 1987A. 

A least squares fit of the tail (from +140 d to +465 
d) yields an e-folding time of ~ 119 days. Comparing the 
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Figure 5. The UBVRI bolomctric light curve of SN 1999em together with that of SN 1987A for comparison. The dotted line corresponds 
to the best fit of SN 1999em late data by SN 1987A tail. The window on the right shows the resulting tail scaled up by an amount of 
0.19 dex (provided by integrated J,H,K fluxes) along with the 56 Co decay for 0.02 of 5B Ni. The reported error bars are la. 



bolometric luminosities of the tail of SN 1999em with SN 
1987A we obtain a best fit of SN 1999em data to the tail of 
SN 1987A for L(99em)/£(87A) = 0.30. With the 56 Ni mass 
0.075 M Q in SN 1987A we thus derive for SN 1999em the 
amount of ejected 56 Ni ~ 0.022 M Q . The bolometric lumi- 
nosity should be corrected for the contribution of the flux 
in JHK bands. According to our photometric observation on 
day 440 the contribution of the JHKs flux is about 0.19 dex. 
Then assuming a similar constant percentage contribution 
by the infrared flux during the tail phase (which is not pre- 
cise), we scale up the derived UBVRI light curve by 0.19 dex. 
The resulting late curve is shown in the window in Figure 5 
together with the best theoretical fit of the radioactive decay 
energy input assuming that the envelope is optically thick 
to 7-rays (L = L x(M m / M e )e- t/TC ° , with i = 1.32x 10 43 
ergs s" 1 and r Co = 111.23 days; Woosley et al. 1989). The 
best fit is found for 56 Ni mass 0.02 ±0.0013 Mq. Both meth- 
ods thus lead to the consistent value of 56 Ni mass of 0.02 Mq . 
This amount is similar to the one derived for SN 1991G (~ 
0.024 Mq; Blanton et al. 1995) and smaller than the mass 
determined for the typical SNe IIP 1969L and 1988A (~ 
0.07 M e ; Turatto et al. 1993). 



3 SPECTROSCOPIC EVOLUTION 

The journal of the spectroscopic observations in the opti- 
cal as well as in the infrared are displayed in Table 2. For 



each spectrum the date (col.l), phase (col.2), range (col. 3) 
and the instrument used (col. 4) are presented. The spec- 
tra have been wavelength calibrated using comparison arc- 
spectra of He-Ar or He-Ne lamps, and flux calibrated using 
observations of spectrophotometric standard stars selected 
from Stone & Baldwin 1983, Baldwin & Stone 1984; Hamuy 
et al. 1992 and Hamuy et al. 1994. 



3.1 Optical spectra 

We present the spectral evolution of SN 1999em which ex- 
tends from ~ +9 d to ~ +642 d since explosion. The redshift 
C2=717 km s _1 of the parent galaxy NGC 1637 has been re- 
moved in all spectra . 

In Figure 6A we show the evolution during the pho- 
tospheric phase. The earliest spectra display a blue contin- 
uum indicating temperatures exceeding 10 4 K with broad 
P-Cygni profiles of the hydrogen Balmer lines, He I lines, 
and the Na I D and Ca II lines. In the first spectrum the 
blue wing of Ha absorption indicate expansion velocities up 
to ~ 16000 km s" 1 . 

Those facts, broad P-Cygni profiles of hydrogen lines, 
high expansion velocities and high temperatures, imply that 
the supernova is a type II SN with a substantial hydrogen 
envelope discovered close to the explosion time. Note that 
the He I 5876 A with a P-Cygni profile is visible for the 
3 and 4 Nov spectra (the same behaviour was seen in SN 
1987A). In fact, in early time spectra Baron et al. (2000) 
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Figure 6. Spectral evolution of SN 1999em during the photospheric phase(A), and during the nebular phase(B). All the spectra have 
been corrected for the recession velocity of the host galaxy NGC 1637 and vertically displaced by an arbitrary constants for clarity. Each 
spectrum is labelled by the corresponding days after explosion. The strongest features are also marked as absorption in the photospheric 
phase(A) and at the rest frame in the nebular phase(B). 



have found evidence for helium enhanced by at least a factor 
of 2 over the solar value as well as a possible nitrogen en- 
hancement. We do not find clear evidence for nitrogen lines 
in our first spectra. As to the He I line, the interpretation 
of its strength in terms of He overabundance is premature, 
since ignored freeze-out effects could lead to the enhanced 
He excitation compared to the steady-state model (Utrobin 
& Chugai 2002). 

At the middle of the plateau phase the photosphere 
cools down to 5000 — 6000 K; at this epoch the spectra (Fig. 
6A) are dominated by strong lines of H, Ca II, Fe II, Na 
I with well developed P-Cygni profiles which become nar- 
rower as the velocity at the photosphere decreases. We note 
also the appearance of absorption lines at 5476 A and 5608 
A which we attribute to lines of Sc II 5526 A and 5658 A, 
respectively although secondary evidence from other lines of 
these ions is not obvious. As the supernova ages, the inter- 
nal energy in the form of the trapped radiation is exhausted, 
the luminosity drops signaling the end of the photospheric 
epoch and transition to the nebular phase. The continuum 
formed in the region excited by 56 Co decay becomes faint 
and diluted compared to black-body intensity. As a result 
the contrast between the net line emission and continuum 
increases, which is demonstrated e.g. by Ha and Ca II 8600 
A triplet (last spectrum in Fig. 6A). 

Figure 6B shows the spectra during the nebular phase, 
from ~ +140 d until ~ +642 d after explosion. The su- 



pernova is now in the radioactive tail phase. The first two 
spectra signal the onset of the nebular phase, where rela- 
tively broad P-Cygni profiles of Ha and Na I 5890,5896 A 
were still seen with FWHM of Ha about 3500 km s _1 . About 
five months later, ~ +312 d after explosion, Ha has become 
still narrower with FWHM of ~ 2400 km s" 1 . The feature at 
7300 A is identified with [Ca II] 7392, 7324 A doublet always 
observed in nebular spectra of SNe IIP. What is unusual is 
the notable domination of the [Ca II] 7300 A over Ha in 
the latest spectra. Note also the emergence of the nebular 
emission of the [O I] 6300,6364 A doublet. Its profile evolu- 
tion can be seen in Figure 6B. The 140 d and 166 d spectra 
showed already the presence of the [O I] although the two 
components were not yet resolved. As time progresses the 
[O I] 6300 A component increases relative to the [O I] 6364 
A one, indicating that the lines become optically thin as the 
supernova envelope expands. 

The spectra of SN 1999em at photospheric and nebular 
epochs in many respects are similar to those of SN 1987A 
(Fig. 7). Note, however, the line width and hence velocities in 
SN 1987A are higher e.g. at ~ +40 d: V%7a{H<x) ~ 6900 km 
s _1 , Vggem(Ha) ~ 6100 km s _1 from the P-Cygni absorption 
minima. The Ba II 6142 A absorption line of an s-process 
element, is very prominent for SN 1987A, while its presence 
is only hinted in the spectrum of SN 1999em. The Sc II 
6280 A line, if correctly identified, on the other hand has 
almost the same strength for the two spectra. The 4554 A 
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Figure 7. The spectral evolution of SN 1999em with that of SN 1987A at similar phases. The spectra of SN 1999em have been shifted 
by 3.5 dex (log-scale). The flux scales correspond to SN 1987A. The spectra are corrected for redshift as well for reddening by the same 
amounts as Fig. 3. 



feature, partly due to Ba II 4554 A, is similar for the two 
SNe. This is possibly due to a blend with other lines, in 
particular Fe II 4555 A and Ti II 4549 A. Note also the 
emergence of the [Fe II] 7155 A emission line. Indeed we 
measure the maximum intensities of this line and comparing 
it with the strong [Ca II] 7291,7324 A line we obtain at 
~ 390 d: J max ([Ca II])/l max ([Fe II]) = 6.98 for SN 1999em, 
while it is around 24.68 for SN 1987A. This suggests that the 
[Fe II] lines do not originate from newly synthesized iron. 

In Figure 8 we report the luminosity evolution of some 
nebular lines, namely Ha, [Ca II] 7291,7324 A (upper panel) 
and [O I] 6300,64 A (lower panel) for SN 1999em. The same 
data for SN 1987A, at similar times, are also shown for com- 
parison. The evolution of Ha and [Ca II] 7291,7324 A are 
clearly decreasing in time for both supernovae in a similar 
and parallel way. This reflects the exponential decline at late 
times. On the other hand the derived luminosities for SN 
1999em are lower than for SN 1987A which must be related 



to the nature and diversities in the progenitor stars of the 
two events. It is noteworthy that the ratio of Ha luminosity 
of SN 1999em and SN 1987A between days 300 and 400 is 
pa 0.3, i.e. the same as the ratio of the 56 Ni. This is perhaps 
not surprising since the rate of Ha emission is determined 
by the rate of radioactive decay. However, we know that 
hydrogen photoionization from the second level is involved 
also in producing Ha quanta (Chugai 1987; Xu et al. 1992), 
which, although it also is eventually determined by radioac- 
tive decay, might spoil the scaling of Ha luminosity with 
56 Ni mass especially at the early nebular epoch. Moreover, 
at the early epoch collisional de-excitation probably plays a 
role also. On the other hand at later epoch (> 400 d) the 
escape of 7-rays should be significant, with the escape rate 
dependent on the mass and explosion energy. Therefore the 
epoch 300 — 400 days is just the right phase at which the Ha 
luminosity scales as 56 Ni mass in SNe IIP. This suggests to 
us the idea of using spectrophotometry of SNe IIP around 
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Table 2. Spectroscopic observations of SN 1999em 



Date (UT; dd/mm/yy) 


Phased 


Range (A) 


Instrument 


* Optical observations : 








03/11/99 


+9 


3240-9060 


Dl 


04/11/99 


+10 


3160-9050 


Dl 


05/12/99 


+41 


3340-10190 


EI 


18/12/99 


+54 


3370-10200 


EI 


29/12/99 


+65 


3340-9050 


Dl 


31/01/00 


+97 


3530-7520 


Al 


1 fi /ri9 /on 


+113 




Dl 


13/03/00 


+140 


3340-9050 


Dl 


08/04/00 


+166 


3380-10300 


EI 


01/09/00 


+312 


3350-10260 


EI 


no /no /fin 

Uo/Uy/UU 


-holy 


Q/1 1 fl 777(1 

o41U- t 1 IV) 


Al 


19/11/00 


+391 


3500-9110 


DI 


01/02/01 


+465 


3360-10260 


EI 


18/03/01 


+510 


3560-9240 


DI 


27/07/01 


+642 


4200-9940 


Ell 


* IR observations : 








14/11/99 


+20 


9420-25040 


EIII 


11/01/01 


+444 


9430-16540 


EIII 



f since the explosion time: ~ 24.5/10/1999 
Al = Asiago 1.82m+AFOSC; DI = Danish 1.54m+DFOSC; 
EI = ESO 3.6m+EFOSC2; EII = ESO-VLT-U1+FORS1 ; 
EIII = ESO-NTT+SOFI 



day 300 in the Ha band with low resolution (say 1000 km 
s _1 ) to measure 56 Ni mass based upon comparison of Ha 
luminosity with the template SN 1987A, without the need 
to reconstruct the bolometric light curve, which would be 
difficult for distant SNe IIP. 

The behaviour of the luminosity of [O I] doublet in both 
supernovae is similar with a "plateau-like" maximum from 
~ 166 d until ~ 391 d for SN 1999em followed by a clear 
steep decline later on, as for SN 1987A. We measured this 
decline rate for SN 1999em to be ~ 0.0166 mag day - . 

Yet the luminosity of [O I] doublet in SN 1999em is a 
factor of 15 lower between day ~ 300 and ~ 500 which in- 
dicates a lower amount of oxygen compared to that of SN 
1987A. To get a rough idea of the ratio of oxygen mass 
in both supernovae one notes that the luminosity of [O I] 
doublet at the epoch of ~ 1 year is powered by the 7-ray 
deposition and by ultraviolet emission arising from the de- 
position of 7-rays in oxygen-poor material. Generally, one 
may write the [O I] doublet luminosity as: 



L(6300)=^^L( 56 Co), 



(1) 



where Mo is the mass of oxygen, M ex is the "excited" mass 
in which the bulk of radioactive energy is deposited, and r\ 
is the efficiency of transformation of the energy deposited in 
oxygen into the [O I] doublet radiation. Given 15 times lower 
luminosity of [O I] doublet in SN 1999em (before dust forma- 
tion) and factor 3.4 lower 56 Ni mass in SN 1999em (section 
2.4), and assuming that in both supernovae r\ and M ex are 
similar, we derive the rough estimate that the mass of oxy- 
gen in SN 1999em is a factor 4.4 lower than in SN 1987A. 
The oxygen mass in SN 1987A according to different de- 
terminations (Fransson & al. 1993; Chugai 1994) is in the 
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Figure 8. The temporal evolution of the luminosity of some neb- 
ular lines for SN 1999em with that of SN 1987A at similar times. 
Upper panel: the [Ca II]7291,7324 A and Ha evolution. Lower 
panel: the [O I] 6300,64 A doublet evolution. Filled symbols cor- 
respond to SN 1999em while open symbols represent SN 1987A 
data. 



range 1.5 — 2 Mq, which translates into ~ 0.3 — 0.4 M§ of 
oxygen in SN 1999em. Given the nucleosynthesis computa- 
tions (Woosley & Weaver 1995) this mass corresponds to the 
main-sequence stellar mass of 13 — 14 M©. 



3.2 Infrared spectra 

The infrared spectra are shown in Figure 9 where some line 
identifications are given. The first spectrum was obtained 
during the photospheric phase (~ +20 d) while the second 
one corresponds to ~ +444 d after explosion during the neb- 
ular phase. The earlier spectrum is dominated by a strong 
continuum with clear evidence of the Paschen series of hy- 
drogen, namely P/3, P7 and P<5 displaying P-Cygni profiles 
as did the optical Balmer lines. The structure of the P7 pro- 
file results from blending with He I 10830 A. The position of 
the peak emission and the absorption minimum of P7 are 
reported in Figure 10, using a linear interpolation between 
the respective positions for PS and P/3. Note the difference 
in velocity. 

The velocity of He I 10830 A absorption component is 
in fact similar to that of Ha estimated by interpolation be- 
tween days 10 and 41. The presence of He I 10830 A at this 
epoch may be a result of the non-thermal excitation due to 
the presence of small amount of 56 Ni in the outer layers of 
the ejecta. Alternatively, excitation can result from the re- 
combination of frozen ionization, the mechanism suggested 
by Utrobin & Chugai (2002). The second spectrum is domi- 
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Figure 9. The infrared spectra on day 20 and 444 after correcting 
them from the redshift of the host galaxy. Some line identifications 
for both phases are also reported. 

nated by emission lines as the absorption components of the 
P-Cygni profiles became weaker. 

Combining these IR spectra together with optical spec- 
tra obtained at similar phases we have extended spectra at 
two quite different phases. Figure 11 displays the resulting 
spectra, after being corrected for redshift of the host galaxy 
as well as for reddening, together with three black-body fits 
to the photospheric phase spectrum corresponding to differ- 
ent temperatures, yielding a temperature of the order 5500 
K as the best fit. Both combinations required no arbitrary 
shifts between optical and IR sections. 

3.3 Ha evolution, "Bochum event" and 56 Ni 
asymmetry 

The evolution of the Ha profile from the early photospheric 
epoch through the developed nebular stage is in many re- 
spects similar to that of SN 1987A: early photospheric 
epoch with normal P-Cygni lines expected in the spherically- 
symmetric case, development in the Ha line of a fine struc- 
ture resembling "Bochum-event" in SN 1987A at the late 
photospheric epoch, and the development of fine structure 
and red-shift in the emission maximum of Ha at the nebular 
epoch (Fig. 12). 

The first two spectra show an undisturbed P-Cygni pro- 
file of Ha with a clear blueshift of the peak emission, in 
contrast with the standard picture of line formation in an 
expanding atmosphere that generally requires an unshifted 
emission maximum. Similar behaviour is also seen in spec- 
tra of other SNe II, namely SN IIP 1988A (Turatto et al. 
1993), SN IIL 1990K (Cappellaro et al. 1995) and SN 1987A 




-2x10" -1x10" 10* 



Velocity (km/s) 

Figure 10. He I 10830 A on day 20. The line is blended with P7 
which strongly affects the He I emission component, but does not 
affect the absorption component of 10830 A line whose absorption 
minimum is shown by a vertical tick at a velocity of S3 11000 km 
s _1 . The expected positions of absorption and emission compo- 
nents of P7 are indicated with the "U" shape segment. Note that 
P7 and He I 10830 A velocities are different. The vertical short- 
dashed line corresponds to the He I 10830 A rest wavelength. 

for which this behaviour has been a subject of theoretical 
studies and has been explained as being due to the reflec- 
tion of photons by the photosphere (Chugai 1988; Jeffery 
& Branch 1990). There is a similar blueshift in the emis- 
sion peak of H/3, clearer for the first spectra. Figure 13(left) 
displays the amount of the blueshift seen in our spectra. 
As the supernova evolves, the blueshift decreases, so that 
by ~ 65 d it approaches zero, more clearly for Ha than for 
H/9 where line blending is a complication. The blueshifts for 
the Paschen hydrogen lines from our first infrared spectrum 
are also shown to be in good agreement with Ha, indicating 
similar properties of the photosphere in the optical as well 
as in the infrared. 

In Figure 13 (right) we display the evolution in time of 
the Ha velocity of SN 1999em together with that of SN 
1987A. Velocities are derived from the minima of the P- 
Cygni absorption components, and in SN 1999em are lower 
than in SN 1987A at similar phases. The nature and the 
similarity of the evolution confirm the initial rapid decline in 
velocity corresponding to the rapid cooling period, followed 
by a slowly decreasing velocity period as expected during 
the plateau phase. The expansion velocities at ~ 97 d and 
~ 113 d and also to a lesser extent at ~ 65 d deviate from 
an otherwise smooth trend. This is due to the emergence 
of some structures in the absorption components that make 
the determination of the minimum absorption complicated 
and doubtful. The evidence of some structures in SN 1999em 
spectra has been noted first by Leonard et al. (2000). 

Around ~ 65 d we note some flattening of the absorp- 
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Figure 11. The resulted extended spectra from the best combi- 
nation of infrared and optical spectra. Both spectra are corrected 
for redshif. The earliest OP(+41 d)+IR(20 d) spectrum, shown 
on the top, is also corrected for reddening. Black-body curves at 3 
different temperatures are plotted with the best fit corresponding 
to a temperature of the order ~ 5500 K . 

tion component with a hint of a weak peak near 6496 A 
in the transition between the minimum and maximum flux. 
This subsequently becomes more pronounced. A less pro- 
nounced feature develops on the red side of the emission 
component at ~ 97 d near 6675 A. All these features show 
a tendency to smaller velocities (i.e. they move toward the 
rest wavelength) as time increases, reaching their maximum 
prominence around +113 d. Similar but not so clear be- 
haviour is evident for other lines in SN 1999em, in particular 
Na I D. Figure 14 illustrates the structure and the presence 
of these "bumps" on the lines at day 113. On the Ha pro- 
file we mark three detectable structures: 1 blue and 2 red 
bumps. In the lower plot we show the Na I D profile which 
reveals clearly the presence of two red bumps at velocities 
very close to ri and 12 seen in Ha. Moreover we note that 
the bi feature velocity is in agreement with the photospheric 
velocities as derived from weak metal lines (Fe II 5018 A and 
5169 A) at corresponding time. For the +113 d spectrum : 
Vp e ii ~ 2100 km s _1 while the corresponding velocity from 
the peak of the bi feature is of the order Vbi ~ 2400 km 
s" 1 . 

This complicated structure of the Ha profile was also 
noted for SN 1987A (Ha fine structure and Bochum event; 
Hanuschick 1988; Phillips et al. 1989; Sartori et al. 1990) 
and in SN 1988 A (Turatto et al. 1993). It occurred earlier in 
SN 1987A and has been ascribed to clumping and mixing of 
radioactive material into the outer envelope with an axially 
symmetric geometry, two-sided jet, by Lucy (1988). This 
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Figure 12. The evolution of Ha spectral region. The spectra 
are corrected for the host galaxy redshift and scaled arbitrary in 
intensity, some spectra are also multiplied by a shown factors for 
clarity. The rest wavelengths and days since explosion are also 
reported . 

possibility is illustrated by Figure 15(left) which shows the 
spectrum on day 97 with the model Ha computed assuming 
a spherical distribution of optical depth and net emission 
source functions. The residual shows excess in the blue and 
red, which conceivably reflects the overexcitation produced 
by 56 Ni jets. It is, of course, not a unique representation of 
the spherical model, but this one at least demonstrates this 
possibility. An alternative view on the fine structure in Ha 
of SN 1987A suggests that the blue bump, on the contrary, 
reflects underexcitation of hydrogen at some velocity in the 
atmosphere related to the deep recombination following the 
cooling wave (Chugai 1991a). Recently this conjecture was 
confirmed by the non-steady state calculations of the hy- 
drogen recombination in the atmosphere of SN 1987A with 
the inclusion of reactions with H~ and hydrogen molecules 
(Utrobin & Chugai 2002). This modelling also confirmed 
the dominant role of the ionization freeze-out effect in the 
population of the second hydrogen level and the Ha forma- 
tion at the photospheric epoch (Chugai 1991b). According 
to the latest computations the non-monotonic behaviour of 
the excitation in the atmosphere of SN 1987A (and SNe IIP 
in general) would be a natural outcome of the freeze-out 
effects, combined with the neutralization of H + in the ion- 
molecular reactions. Of course, only direct similar modelling 
for SN 1999em will give a conclusive answer concerning the 
plausibility of this mechanism for the blue bump. Figure 
15(right) demonstrates, how the Sobolev optical depth in 
Ha would appear in order to reproduce the appropriate blue 
bump. Note, although the blue structure is fitted quite well, 
some excess in the red remains unaccounted for in the model 
and therefore indicates a real asymmetry of the hydrogen 
excitation. 
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Figure 13. Left panel: the amount of the blueshift seen in the emission peaks of Ha and H/3 during the first phases. Line blending is a 
complication for H/3. Right panel: the evolution of the expansion velocities, derived from the absorption minima in Ha profiles, for SN 
1999em together with that of SN 1987A for comparison. For both panels the corresponding velocities of P/3 and P<5 from the ~ +20 d 
infrared spectrum are also shown. 
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Figure 14. The Ha and Na I D line profiles at ~ +113 d after 
explosion. The rest wavelengths are also shown. The short bars 
in Na I D profile correspond to the derived velocities from the 
features seen in the Ha profile (see discussion in the text). 

In fact soon after, on day 113 the red excess became 
even more evident. It increases with time and on day 444 it 
is also detected in He I 10830 A. On days 140 and 166 the Ha 
maximum shows significant deviation from a round-shaped 
form. The complicated structure with the dominant red peak 



apparently indicates the asymmetric distribution of the line- 
emitting gas in the central region v < 1500 km s^ 1 . The 
comparison of He I 10830 A with Ha at the similar epoch 
(Figure 16) shows that the peak is redshifted by the same 
value as in Ha. We believe that this red bump at ~ +400 
km s _1 is caused by the 56 Ni asymmetry. The He I line is 
narrower (in the red P7 contributes) . The blue width at half 
maximum (BWHM) is 2000 km s" 1 for Ha and only 1000 
km s _1 for He I 10830 A. The He I 10830 A line traces non- 
thermal excitation more closely than Ha, since for hydrogen 
additional ionization from the second level is more important 
than for He I. The width and position of the He I 10830 A 
line thus indicates that the bulk of 56 Ni is distributed inside 
a sphere with velocity v < 1500 km s _1 and the 56 Ni zone 
is shifted to the far hemisphere by roughly 400 km s _1 . 



4 EVIDENCE FOR DUST FORMATION 
AROUND DAY 500 

The temporal evolution of the [O I] profile can be an effective 
indicator of dust formation in the form of an observable 
blueshift in the emission peaks of the components (Danziger 
et al. 1991). It should be noted, that apart from the late 
time blueshift related to the dust formation, SN 1987A also 
showed an "early blueshift" in [O I] 6300 A between day 200 
and 300 that disapeared at the later epoch. We believe, the 
blueshift in [O I] 6300 A observed in SN 1999em on day 312 
and 391 and which vanished on day 465 is of the same origin 
as the early blueshift observed in SN 1987A and also that 
reported for SN 1988H (Turatto et al. 1993). Namely, it has 
nothing to do with the dust and probably is the result of 
the superposition of the blend of Fe II lines (multiplet 74). 
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Figure 15. Left panel: the Ha on day 97. The overplotted (thick line) is a profile for a spherical model. The residual (bottom panel) 
shows two peaks, which can be interpreted as evidence for the presence of bipolar excitation regions in the atmosphere. Right panel: the 
H« on day 97. The overplotted (thick line) is the profile for the spherically symmetric model with the non-monotonic behaviour of the 
optical depth shown in window (up right corner). The residual shows the excess in the red, which is indicative of overexcitation in the 
far hemisphere. 



This latter is hinted at by the bump at 6250 A attached to 
the blue wing of [O I] 6300 A. 

However between day 465 and 510 the [O I] 6300,6364 
A line and Ha demonstrate a pronounced evolution (Figure 
17). The transformation in both profiles can be described 
as a flattening of the line accompanied by skewing towards 
the blue. The quality of the 642 d spectrum is not sufficient 
to confirm this effect, although the spectrum does not con- 
tradict the conclusion from the previous two spectra. The 
dramatic change of the line profile at the time scale of the 
order of 0.1 of the expansion time (or even faster) and a 
proximity of the age to the epoch of the dust formation in 
SN 1987A around day ~ 526 (Lucy et al. 1989) suggest that 
dust formation is the likely cause of this profile transforma- 
tion. Note that the flattening of [O I] 6300 A line profile 
in fact was not observed in SN 1987A at the dust forma- 
tion epoch. Instead the blueing of oxygen line in SN 1987A 
preserved the round-topped profile of the [O I] 6300 A line 
which is consistent with the dust formation in the sphere 
with the expansion velocity of 1800 km s _1 and dust optical 
depth of the order of unity (Lucy et al. 1989). 

Nevertheless, the flat-topped profile is expected, if the 
optical depth of the opaque core is very high. To emphasize 
the point we present the toy model (Figure 18), which sug- 
gests some smooth distribution of emissivity in the [O I] dou- 
blet lines (with the line ratio of 0.4, similar to SN 1987A at 
the same epoch) and the embedded opaque core. The model 
doublet in the absence of any absorption roughly represents 
the situation on day 465. Switching on the dust absorption 
in the dusty core with the boundary velocity of 800 km s _1 
produces the effect, which markedly depends on the optical 
depth. For r = 2 the line shape and blue shift recalls what 
was seen in SN 1987A at the dust formation epoch after day 



530, where the estimated optical depth of the core is r d ~ 1 
(Lucy et al. 1989). We see that the flat-topped [O I] doublet 
profile observed in SN 1999em on day 510 requires an optical 
depth of the dusty sphere Td 3> 10, since in case of Td = 10 
the plateau has a noticeable inclination, which is absent in 
the case of r d = 50. This illustration indicates that the dust 
formed between days 465 and 510 and the optical depth of 
the dusty core was enormous, rd 3> 10. The blue end of the 
plateau in the [O I] 6300 A line (Figure 17) indicates the 
velocity of the dusty sphere of about t>d ~ 800 ± 50 km s _1 . 

As in SN 1987A, for the Ha line the effect of the dust 
core is less apparent because of the more extended Ha- 
emitting zone and, primarily, because of the redward asym- 
metry of this line related with the 56 Ni asymmetry. Never- 
theless, the blueward skewing in the spectrum on day 510 is 
seen in Ha also. 

Even less apparent is the profile evolution of the [Fe II] 
7155 A line. The absence of the skewing towards the blue in 
this line on day 510 indicates that the unocculted fraction of 
asymmetric red component is strong enough to maintain the 
original redshift of this line seen on day 465. To illustrate 
this explanation we represent the profile of [Fe II] 7155 A as 
a combination of the spherically-symmetric component and 
the component originating in the conic structure with the 
opening angle 2a and the angle between line of sight and the 
cone axis 6 (Figure 19, upper panel). Adding to the model 
the dust core with v d = 800 km s _1 and r d = 50, we checked, 
if our choice for the asymmetrical component is tolerated by 
the 510 day spectrum of [Fe II] 7155 A . The decomposition 
of the line profile into symmetric and asymmetric component 
is, of course, not unique, given the large noise in the 510 d 
line (the usual drawback for the inverse problem). However, 
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Figure 16. Ha and He I 10830 A lines at the nebular epoch. 
Both lines show apparent redshift indicating the asymmetry of 
56 Ni. The He I line blue width at half maximum is factor two 
lower, possibly indicating the actual width of the 56 Ni zone. Note 
that the He I 10830 A maximum is shifted arbitrary to reach Ha 
maximum for comparison. 

we found that for small inclination angle, 9 < 50° , the skew- 
ing towards the blue in the model is unacceptably large. One 
of the acceptable possibilities found by trial and error proce- 
dure is shown in Figure 19(lower panel). In this case 9 = 67° 
and a = 40°. The line emissivity in the model cone is con- 
centrated towards the cone axis and distributed along the 
radius with the broad maximum in the range 1000 — 3000 
km s _1 . In this modelling we followed a somewhat arbitrary 
but sensible requirement that the radial displacement of the 
maximum of the jet emissivity is minimal. The model with 
the dusty core obviously reproduces the major property of 
the profile on day 510, namely, the shift of the maximum 
towards the red. We thus conclude that the dust core im- 
plied by the [O I] doublet transformation is not in conflict 
with the [Fe II] 7155 A line profile provided that the 56 Ni 
jet, responsible for the red asymmetric component, has rel- 
atively large inclination angle 9 > 60°. We do not insist on 
the uniqueness of the model of asymmetric distribution of 
emissivity of [Fe II] 7155 A. Yet the one-sided asymmetry of 
emissivity implied by the profile evolution is strikingly con- 
sistent with the one-sided 56 Ni distribution suggested above 
by Ha and He I 10830 A lines. 

The manifestations of the dust formation in SN 1999em 
are different from those in SN 1987A in at least three re- 
spects. First, while in SN 1999em the [O I] 6300 A line 
revealed the presence of the dust already on day 510, in 
the case of SN 1987A the blue shift in [O I] 6300 A 
emerged only after day 526 (Danziger et al. 1989; Phillips & 
Williams 1991). A somewhat earlier emergence of the dust 
in SN 1999em is probably related to the lower amount of 
56 Ni, which in its turn implies lower temperature at the 



Figure 17. The oxygen doublet [O I] 6300, 6364 A and Ha in 
the latest three epochs. The spectra on day 642 are multiplied by 
factor of six. The two vertical lines correspond to zero velocity 
and the approximate position of the blue edge of plateau in [O I] 
6300 A line. The spectra are not corrected for reddening. 

given stage, so the condensation temperature was attained in 
SN 1999em earlier. Second, the velocity of the dust-forming 
zone is a factor two lower compared to the dusty core in 
SN 1987A (~ 1800 km s _1 ) measured in the spectra around 
day 600 (Lucy et al. 1989). This difference possibly reflects 
both the smaller mass (and therefore velocity) of the metal- 
rich core of ejecta and smaller velocity of 56 Ni bubble, which 
presumably pressurizes the metal-rich zone thus produc- 
ing dense clumps in which dust forms. The third is more 
clear cut: the optical depth of the dusty core in SN 1999em 
r d > 10 while in SN 1987A r d < 1 (Lucy et al. 1991). At 
first glance a factor two smaller radius of the dust core in 
SN 1999em compared to SN 1987A is an explanation for this 
difference. Assuming a comparable amount of dust and the 
given factor of two smaller radius one gets a factor of four 
higher column density. However, this is still a small factor to 
explain at least a factor of 20 difference in the optical depth. 

The explanation of the high optical depth derives from 
the model of the clumpy dust zone in SN 1987A proposed by 
Lucy et al. (1991). To account for the wavelength indepen- 
dent extinction the authors suggested that the dust is locked 
in very opaque clouds, so the effective optical depth of the 
dusty zone is actually the geometrical (occultation) optical 
depth produced by the cloud ensemble (r oc = 3N/4irRa, 
where N is the number of clouds, i?d = v^t is the radius 
of dusty core). Lucy et al. (1991) estimate r oc ~ 0.4 for 
SN 1987A around day 625. To account for the high optical 
depth of the dusty core in SN 1999em, we must admit that 
the occultation optical depth is very high r oc 3> 10, i.e. the 
number of opaque cloudlets is substantially higher than in 
SN 1987A. We estimate the amount of dust required to pro- 
duce the optical depth rd = 50 on day 510 as ~ 10 -4 Mq, 
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Figure 18. The effect of the opaque core on the [O I] doublet. 
Lines of different thickness correspond to different optical depths 
of the core, starting with r = (thinnest line) through t = 2, 
t = 10, and r = 50 with progressively growing line thickness. 
Note, for r = 10 the plateau is still not settled and only for 
r 2> 10 the profile is flat-topped. 



a quite moderate value. In fact, this number should be con- 
sidered as a lower limit. An apparent blueward shift in the 
peak of Mg I] 4571 A is apparent between days 465 and 510. 
Since the shift is similar in velocity to that observed in [O 
I] 6300 A and clearly results from a depletion of flux on the 
red side, with the blue side remaining constant in velocity, it 
is consistent with our model for a very opaque core of dust 
with a large occultation optical depth. 

The opaque dusty core must radiate as a black-body 
with the temperature determined by the absorbed luminos- 
ity and the radius. Modelling the Ha profile transformation 
we estimate that about 30% of the luminosity is intercepted 
by the dusty core on day 510. Given the bolometric lumi- 
nosity at this epoch L ~ 2 x 10 39 erg s _1 the luminosity 
of the dusty core must be then of Ld ~ 6 x 10 38 erg s _1 . 
The black-body temperature of the dusty core on day 510 
was then pb 510 K with the maximum of the infrared spec- 
trum at 6 nm. We are not aware of any IR observations 
of SN 1999em in M band at this epoch. Since vibration- 
rotation bands of CO have been detected in the IR spectra 
of SN 1999em by Spyromilio et al. (2001), one is entitled to 
conclude, as was done with SN 1987A, that the formation 
of molecules was a prelude to formation of dust. Moreover 
we expect dust condensation to cause a deficit in the opti- 
cal radiation with respect to the case with no dust. Figure 
20 displays an expanded version of the B, V and R late 
time light curves together with the least squares fits (the 
last point is excluded in the fit). It is clear that all the light 
curves display a deficit decline for day 510 compared with 
the linear trend shown by fitting the previous data as repre- 
sented by the dashed line. The average deficit is close to 30 
percent. This gives further support for dust formation after 
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Figure 19. The effect of the opaque dusty core on the [Fe II] 7155 
A line. Top panel shows the model line profile model (thick line) 
overplotted on the observed profile on day 465 (dotted line). The 
model is a combination of the spherically-symmetric and asym- 
metric components (thin lines). The asymmetric component is 
modeled by the conic overexcitation zone in the far hemisphere. 
The lower panel shows a similar plot but on day 510 and with 
added opaque core velocity 800 km s — 1 and r = 50 in the model. 
Note that the model flux on day 465 was multiplied by factor 0.5 
to produce the adequate fit on day 510. 
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Figure 20. The late B,V and R light curves. The dashed lines 
represent the least squares fits of data weighted according to error 
estimates and excluding the last point at day 510. A deviation 
is seen in the last point of each colour, which supports the dust 
condensation scenario after day 465 as inferred from spectroscopy. 
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day 465. Such a deficit was also noted by Lucy et al. (1991) 
at the time of dust formation in SN 1987A. 

The fact that a fully opaque dusty core absorbs only 
a relatively small fraction (~ 30%) of the optical radiation 
is apparently related to the compactness of the dusty core 
and the large extent of the optical emitting zone. This is 
indicated by the fact, that the velocity of the dusty sphere is 
notably lower than the FWHM of Ha on day 465 (» 2700 km 
s _1 ). A simple geometrical model to illustrate the situation 
is the dusty sphere with the velocity embedded in the 
homogeneous extended emitting sphere (v c ). It is easy to 
show then that the absorbed energy fraction ~ 0.3 implies 
the ratio Vd/v c ~ 0.6. Given va ~ 800 km s^ 1 , this leads to 
the effective velocity of the emitting sphere of v c ~ 1300 km 

a" 1 . 

We may use a more physical model to directly constrain 
the extent of a 56 Ni zone (vm) from the relative amount of 
absorbed light. Let the SN envelope be a homogeneous freely 
expanding sphere with the outer velocity vo determined by 
the ejecta mass M and kinetic energy E. We assume that 
56 Ni is distributed homogeneously in the range v < vm- We 
may then compute the bolometric luminosity (without in- 
frared dust emission) treating the gamma-ray transfer in the 
absorption approximation with the conventional effective 
absorption coefficient k = 0.03 cm 2 g _1 . The positron depo- 
sition is assumed to be local. Simulations with and without 
dust using vm as a tuning parameter permit us to find the 
extent of the 56 Ni zone from a condition that the dusty core 
with Vd = 800 km s _1 and ra = 50 absorbs 30% of the light 
generated locally with the deposition rate. In the particular 
case of the ejecta mass of M = 12 M , and two values of 
kinetic energy of E = 10 51 erg and E — 5 x 10 50 erg we find 
that the required radius of the 56 Ni zone must be «Ni = 1100 
km s _1 and 1300 km s _1 , respectively. The parameter vari- 
ation probably would not change the results markedly, al- 
though an asymmetry could slightly affect the situation. The 
resulting extent of the 56 Ni zone (~ 1100 - 1300) km s _1 is 
quite sensible and consistent with our previous finding that 
56 Ni lies inside the sphere of ~ 1500 km s _1 . Remarkably, 
the above simple geometrical analysis led to the extent of 
the emission region (1300 km s _1 ) comparable to the ex- 
tent of 56 Ni zone in the dense (low energy) model. This 
reflects a simple truth that a dense model is close to the 
situation of the local deposition, the approximation implic- 
itly assumed in the simple geometrical model. We therefore 
conclude that the small partial blackout of the supernova 
light following the formation of an extremely opaque dusty 
core of SN 1999em stems from the fact that the dusty core 
is compact compared to the extended distribution of 56 Ni. 



5 THE PHOTOSPHERIC TEMPERATURE, 
EXPLOSION TIME AND DISTANCE 

At the early photospheric epoch the photosphere of SNe IIP 
is considered to approximate a black-body with a sharp 
boundary (Branch 1990). We made black-body fits to our 
spectra, after dereddening the fluxes in different bands with 
E(B — V)=0.1. The interstellar extinction laws of Cardelli 
ct al. (1989) were adopted. The derived temperatures are 
shown in Figure 20(upper panel), plotted as a function of 
time since explosion. We have used our Nov 3 spectrum to 



estimate the difference in temperature obtained by using a 
reddening of E(B — V)=0.1 and the smaller reddening used 
by Baron et al. (E(B — V)=0.05). That variation is then ap- 
plied to the Baron et al. estimate for Oct 29 spectrum to 
give a revised higher temperature of T~14300 K. In Figure 
21(left lower panel), we plot the evolution of the difference 
in temperature between the case of no reddening correc- 
tion (E(B - V)=0.0) and the one with E(B - V)=0.1. This 
behaviour shows the effect of the reddening correction on 
the derived temperatures from black-body fitting. At early 
phases, close to the explosion time and when the object is 
very hot, the colours on the Ray leigh- Jeans tail are not sen- 
sitive to temperature and therefore small changes in colour 
produce large changes in temperature and as time progresses 
this difference decreases. Figure 21(left upper panel), on the 
other hand, confirms the strong initial cooling of the pho- 
tosphere as the supernova ages followed by a relatively con- 
stant temperature starting at T~5500 K and dropping to 
T~4200 K. 

We apply the expanding photosphere method (EPM) to our 
data in order to estimate the time of explosion of SN 1999em 
and its distance. More detailed studies concerning the use of 
this method in the case of SN 1999em were made by Hamuy 
et al. (2001) and by Leonard et al. (2002), who used broad 
band colours rather than spectrophotometry, and a larger 
data base. The EPM, also called Baade-Wesselink method, 
is based essentially on the assumptions that the SN IIP pho- 
tosphere is spherical and emits as a dilute black-body, and 
that the expansion proceeds homologously (v = r/t). For 
small redshift (z « 1) the photospheric angular radius is 



R 
D 



/a10°-' 



(2) 



where R — vt is the photospheric radius determined by the 
expansion time t and velocity at the photosphere v (here we 
assume that the presupernova is a point) , D is the distance 
to the SN, f\ is the apparent flux density, A\ is the total red- 
dening coefficient, B\(T) is the Planck function evaluated 
at T, and £a is the dilution parameter to account for the 
fact that the SN does not radiate as a perfect black-body. 

We adopt A v = 0.31, fx at A = 5500A (from obser- 
vations), B\(T) using the temperatures of the black-body 
fitting of the spectra, and the parameter £ according to tem- 
perature and wavelength dependence given by to Hamuy et 
al. (2001). The expansion velocity v was measured from our 
spectra, using the absorption minima of weak lines such H/3 
and Fe II lines at 5018 A & 5169 A, formed in material close 
to the photosphere (Eastman & Kirshner 1989). The veloc- 
ity evolution is shown in Figure 21(right). This leaves two 
unknowns, to and D, to determine. Results are reported in 
Table 3. Figure 22 shows a least squares fitting of our data. 
We found to ~ JD 2451476(±lday) corresponding to 1999 
October 24.5 and a distance of 7.83 ±0.3 Mpc. These derived 
values agree well with results of the more detailed studies 
of the EPM for the case of SN 1999em (Hamuy et al. 2001 
and Leonard et al. 2002). Our use of the same values of £ 
and the reddening as those used by Hamuy et al. links the 
two results. It is the use of different type of photometry 
and independantly determined velocities that strengthens 
the conclusions by both groups that the results are robust. 
Thus the case of SN 1999em provides a check of the con- 
sistency of the EPM since, from observations, it is already 
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Figure 21. Left panel: the evolution in time of the derived temperatures from black-body fitting of the spectra (A), and the variation of 
the change in temperature from E(B — V)=0.0 to E(B — V)=0.1 (B). Right panel: the photospheric velocity and the photospheric radius 
evolution. This later is computed using the temperature and the bolometric luminosity (sections 2 & 5). 



Table 3. The EPM results for SN 1999em. 



Date (JD 2451000+) 


Flux (10 _14 ergs/s/cm 2 /A in V) 


T bb (K) 


V (km/s) 


t 


6 (10 13 cm/Mpc) 


481.2 


1.101 


14312 


10990 


0.45 


6.034 


485.863 


1.398 


12802 


9440 


0.424 


8.278 


486.767 


1.385 


11217 


9072 


0.42 


9.685 


517.795 


1.228 


5510 


4050 


0.65 


20.707 


530.755 


1.184 


5080 


3280 


0.75 


21.578 


540.7 


1.152 


5040 


3120 


0.8 


20.361 


575.42 


0.874 


4520 


2150 


0.88 


21.772 



The values of £ are taken from Hamuy et al. (2001). 



known that the supernova exploded between Oct 20.45 and 
the discovery date Oct 29(IAUC 7294). In addition the de- 
rived distance is in excellent agreement with that derived 
from a study of the stellar component of NGC 1637 by Sohn 
and Davidge (1998). 

We compute as well the photospheric radius using the 
luminosity and the temperature for a given time. The re- 
sults are shown in Figure 21 (right), together with the pho- 
tospheric velocity evolution. 



6 PROGENITOR STAR PROPERTIES 

The detailed observations of SN 1999em provide us with 
three important phenomenological parameters: the plateau 
duration (t p ), absolute V magnitude at the plateau (Mv) 
and the velocity at the photosphere measured from weak ab- 
sorption lines (Vph). These three parameters permit us to re- 



cover the explosion energy (E), ejecta mass (M) and presu- 
pernova radius (R) using the light curve models as originally 
proposed by Litvinova & Nadyozhin (1985). In this paper 
useful expressions are given, which are used to estimate SN 
parameters. Similar, but somewhat different expressions are 
given by Popov (1993) based on the analytical model, which 
will also be used here. The input parameters of the model 
are the duration of the plateau phase (t p ), a representative 
photospheric velocity at the beginning of the plateau phase 
(V'ph) and its absolute V magnitude {My). From observa- 
tions of SN 1999em we find V p h ~ 4050 km s^ 1 using weak 
lines in the spectrum on day 41, Mv ~ —15.76 at the same 
epoch and the duration of the plateau t p ~ 80 days. Accord- 
ing to the relations from Popov (1993) we derive the follow- 
ing parameters for SN 1999em: ejecta mass M « 11 Mq, 
explosion energy E » 1.1 x 10 51 erg, presupernova radius 
R » 120 Rq. Applying the expressions derived by Litvinova 
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Table 4. Decline rates for the second plateau feature. 



Supernova 


B band 


V band 


R band 


I band 




(mag) 


(mag) 


(mag) 


(mag) 


1997D 


~ 0.33 in 221 d 


~ 0.69 in 89 d 


~ 0.57 in 79 d 


~ 0.53 in 51 d 


1991G 




~ 0.3 in 30 d 


~ 0.2 in 40 d 


~ 0.15 in 40 d 


1999em 


~ 0.22 in 58 d 


~ 0.5 in 58 d 


~ 0.45 in 58 d 


~ 0.54 in 58 d 



The B-light curve of SN 1991G is not available in the work 
of Blanton et al. (1995). 
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Figure 22. The derived EPM results from our data. Also 
shown the least squares fit and the estimated explosion time and 
the distance to the SN. The discovery time corresponds to JD 
2451480.94. 



& Nadyozhin (1985) using a grid of hydrodynamic models 
for type II SNe, we find: ejecta mass M « 9.8 Mq, explosion 
energy E » 0.5 x 10 51 erg, presupernova radius R ~ 180 Rq. 
Both estimates are then consistent with a claim of a pre- 
supernova radius of ~ 120 — 150 -Rq, an ejecta mass of 
10 — 11 Mq and an explosion energy ~ (0.5 — 1) x 10 
erg. With approximately 1.5 Mq enclosed in the neutron 
star and 0.5 — 1 Mq lost by the wind, we thus obtain the 
12 — 14 Mq as a possible range for the main sequence mass of 
the progenitor. This is consistent with the values based upon 
oxygen mass estimate (section 3.1), and as well with the up- 
per limit on the main-sequence mass using pre-supernova 
field images (Smartt et al. 2001). 

The parameters derived using both sets of analytical 
expressions and the radius in particular (~ 120 — 150 Rq) 
suggest that the progenitor was a Go - G5 supergiant (Cox 
2000). If it were at the detection limit of my = 23.2 sug- 
gested by Smartt et al. (2001) it would have had My = 
—6.58 and M& j = —6.42 and an effective temperature of 
the order 5570 K. Using these parameters and the theoreti- 
cal tracks of Maeder & Meynet (2001) we conclude that the 



most probable mass would be 12±1 Mq. This conclusion is 
not significantly affected by uncertainties in rotational veloc- 
ities and metallicity. These estimates provide one with the 
main parameters of the SN progenitor as a starting point 
for more detailed computations. 

An interesting feature of the bolometric evolution of SN 
1999em is the short duration of the peak before the fall to 
the plateau phase. In fact it has been shown that this early 
peak is very sensitive to the presupernova mass-loss history. 
Indeed Falk & Arnett (1977) have demonstrated that a rapid 
rise and short duration peak (~ 1-2 days) is consistent with 
the absence of an extended dense circumstellar shell and 
hence a low mass loss rate immediately prior to the explo- 
sion. This is also in agreement with Chandra x-ray and radio 
observations of SN 1999em (Pooley et al. 2001). 



7 SUMMARY AND CONCLUSIONS 

We have presented photometric and spectroscopic data for 
SN 1999em, from ~ 9 d until ~ 642 d after the explosion. 
The shape of the light curve (t p ~ 80 days) as well as spectral 
features show that it is a type IIP supernova. The problem 
of reddening has been discussed and we conclude that an 
optimum choice is close to that determined by Hamuy et 
al. (2001). The analysis of late phase photometry, up to ~ 
510 d, shows that the exponential tail decay rate is close 
to the one of the radioactive decay 56 Co to 56 Fe, indicating 
that this is the main source of energy powering the light 
curve. 

A photometric comparison of SN 1999em with SN 
1987A, especially in the later phases, provides constraints 
on the radioactive 56 Ni mass. We have constructed the "UB- 
VRI" bolometric light curve, and comparing it with that of 
SN 1987A we obtain an estimate of the amount of 56 Ni pro- 
duced by the explosion of 0.02 Mq, a smaller value than 
that derived for typical type IIP SNe such as SN 1969L and 
SN 1988A which have M( 56 Ni) ~ 0.07 Mq. 

We have noticed some flattening in the light curves, just 
after the steep decline from the plateau phase, and clearer for 
the blue bands. This behaviour is also seen in the light curves 
of two other objects, namely the peculiar SN IIP 1997D 
and the SN IIP 1991G. All three events have some com- 
mon features, being type IIP and all having a lower ejected 
56 Ni mass: a similar amount for both SN 1999em and SN 
1991G,(~ 0.02 M©) while SN 1997D ejected an even lower 
mass (~ 0.002 Mq). In addition, the duration of this "second 
plateau" on the tail seems greater for SN 1997D than for SNe 
1991G and 1999em. The decline rates and the observed du- 
ration of the flattening period are reported in Table 4. These 
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Supernova 


L p 


E 


R 


M( 56 Ni) 




Source 




(days) 


(10 51 ergs) 


(Rq) 


(Mq) 


(Me) 




1999em 


~ 80 


~ 0.5 - 1 


120 - 150 


~ 0.022 


12 - 14 


A 


1987A 


~ 40 


~ 1.3 


~ 40 


0.075 


~ 20 


13 


1997D 


~ 50 


0.1 


85 


0.002 


8 - 12 


C 


1993J 




1.6 




0.078 


12 - 16 


D 


1969L 


~ 90 


~ 1.7 


~ 220 


~ 0.07 


~ 20 


E 



A: The present work D: Utrobin 1996; Chugai & Utrobin 2000 

B: Danziger et al. 1988; Wooslcy et al. 1989 E: Arnett 1996; Sollcrman ct al. 1998 
C: Chugai & Utrobin 2000; (see Turatto et al. 1998 for alternative model) 



measurements suggest that the "second plateau" feature is 
a common feature for this low 56 Ni mass type IIP super- 
novae, and that its duration is correlated with the amount 
of ejected Ni. Further support for this possible correlation 
comes from the case of SN IIP 1999eu (Pastorello et al. in 
preparation) which shows a very clear second plateau feature 
of duration ~ 200 d (clear in V and B bands). Moreover it 
seems that this SN has a very low ejected 56 Ni mass similar 
to SN 1997D. Note that the prototype type IIP supernovae 
SN 1969L and 1988A do not show clear evidence of this be- 
haviour. On the other hand SN 1969L and SN 1988A are 
known to produce an amount of 56 Ni similar to SN 1987A 
~ 0.07 Mq. Improved statistical samples and better sam- 
pled light curves are required in order to confirm or rule out 
this behaviour. Radiation diffusion effects, still important at 
the beginning of the radioactive tail, are a possible cause of 
this behaviour. 

SN 1999em provided a test of the validity of the EPM 
since we have observational constraints on the explosion 
time and distance. In fact analysing the spectrophotome- 
try, we derive an explosion time consistent with these con- 
straints and in good agreement with what was found from 
more detailed broad band photometric studies of SN 1999cm 
(Hamuy et al. 2001; Leonard et al. 2002). 

We analysed the phenomenon of the fine structure of 
Ha at the photospheric epoch, which was reminiscent of the 
"Bochum event" in SN 1987A. Two possible explanations, a 
bi-polar jet proposed by Lucy (1988) for SN 1987A and un- 
derexcitation of hydrogen combined with 56 Ni asymmetry, 
are discussed. This analysis does not permit one to discrim- 
inate between those models. Yet the one-sided 56 Ni ejection 
seems to find support in the small red shift of the He I 10830 
A profile and the larger red shift of the Ha line profile at the 
nebular epoch. We note that the He I 10830 A line should 
be more sensitive to non-thermal excitation resulting from 
7-ray deposition. These lines indicate that the 56 Ni distribu- 
tion could be imagined to be a filled sphere with a velocity 
of ~ 1500 km s _1 shifted to the far hemisphere by 400 km 
s _1 . A somewhat surprising coincidence is that in SN 1987A 
the 56 Ni distribution also shows one-sidedness with a shift 
to the far hemisphere. 

Analysing the [O I] 6300,6364 A line profile evolution 
we found a rapid change between days 465 and 510, which 
we interpret as an effect of the dust formation during this 
interval. Other support for dust formation comes from the 
deficit seen in optical radiation measured by late time pho- 
tometry. It is the second SN IIP (after SN 1987A) where 



convincing evidence of dust formation exists. In SN 1987A 
and SN 1999em we detected the rapid (At/t ~ 0.1) trans- 
formation of profiles during the nebular epoch. The dust 
phenomenon in SN 1999em has some distinctive character- 
istics compared to SN 1987A. The dust condensation hap- 
pened earlier (between days 465 and 510) than in SN 1987A 
(after day 526), which is probably explained, by the lower 
56 Ni mass and, accordingly, lower temperature. The dust 
resides in the core with a velocity of w 800 km s~\ much 
lower than in SN 1987A. This greater confinement of dust 
in SN 1999em possibly results from the lower velocity cre- 
ating a more confined metal-rich region where the condition 
for dust formation prevails. Another remarkable difference 
is the very large optical depth of the dusty zone (r >• 10) 
compared to SN 1987A (r ~ 0.5, Lucy et al. 1991). We in- 
terpret this difference as an indication that in SN 1999em 
the dust is distributed more homogeneously (or the num- 
ber of the opaque dusty clumps is notably greater) than in 
SN 1987A. 

These facts show clearly the importance of studying 
more samples of SNe IIP (SNe II in general), because of 
the diversity they provide in manifesting the same event 
(i.e. dust condensation) and thus the opportunity of under- 
standing the physics behind such events. 

We used relations by Litvinova & Nadyozhin (1985) 
and Popov (1993) to find SN parameters from observa- 
tional characteristics. The estimated mass of the progenitor 
M ms « 12 — 14 Mq and presupernova radius (120 — 150 R@) 
are just consistent with the failure to detect a progenitor 
star with imaging of the pre-SN field. Moreover, our derived 
progenitor mass agrees with our finding that the oxygen 
mass in SN 1999em is about four times smaller compared 
to SN 1987A. 

There is growing observational evidence of a significant 
decrease in the mass of Fe produced as the progenitor mass 
for type II SNe decreases when we include derivations for 
other objects. It is therefore encouraging to note that this 
behaviour is what is required to explain abundance patterns 
such as [O/Fe] and [Mg/Fe] in metal-poor halo stars modeled 
by Argast et al. (2001). 

The correlation of progenitor mass of SN IIP and core 
collapse SN in general with other parameters is vital for test- 
ing explosion models and the theory of stellar evolution. In 
Table 5, we present parameters of SN 1999em along with 
other SNe II. Note that the parameters reported in the ta- 
ble, except for SN 1999em, are not obtained directly by using 
the analytical models of Litvinova & Nadozhin (1985) and 
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Popov (1993). They are however the most reliably deter- 
mined through modelling of observations. 

Table 5 demonstrates that for type II SNe there begins 
to emerge a monotonic relation between progenitor mass on 
the one hand and both explosion energy and 56 Ni mass on 
the other. Although the uncertainties in these parameters 
remain, of necessity, large, it seems that in the progenitor 
mass range 10 — 13 Mq there is also a steep decrease in 
both energy and 56 Ni mass. This suggests that these latter 
2 parameters are not independent but physically linked. Re- 
sults for other core collapse objects (type Ib,c) with higher 
progenitor masses tend to support this correlation even if 
the current scatter is unavoidably large. Future observations 
supported by modelling will surely elucidate this conclusion. 
It is clear from our analysis that high S/N spectra of SNe 
at late phases are an invaluable tool for understanding type 
II supernovae. 

The fact that SN 1999em probably resulted from the 
explosion of a Go — G5 supergiant also indicates, as did SN 
1987A, that the evolutionary stage at which massive stars 
can explode is not yet well delineated. 
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